A new type of limit theorems for the one-dimensional 

quantum random walk 



By Norio KONNO 



Abstract. In this paper we consider the one-dimensional quantum random walk X% at time n 
starting from initial qubit state ip determined by 2 x 2 unitary matrix U. We give a combinatorial 
expression for the characteristic function of X%. The expression clarifies the dependence of it on 
components of unitary matrix U and initial qubit state <p. As a consequence, we present a new 
type of limit theorems for the quantum random walk. In contrast with the de Moivre-Laplace limit 
theorem, our symmetric case implies that Xn/n converges weakly to a limit as n — > oo, where 
Z 1 ? has a density l/vr(l - x 2 )y/\ - 2x 2 for x G (-l/\/2, 1/V2). Moreover we discuss some known 
simulation results based on our limit theorems. 



1 Introduction. 

The classical random walk on the line is the motion of a particle which inhabits the set of 
integers. The particle moves at each time step either one unit to the left with probability 
p or one unit to the right with probability q — 1 — p. The directions of different steps are 
independent of each other. This classical random walk is often called the simple random 
walk. For general random walks on a countable space, there is a beautiful theory (see Spitzer 
|12j). In the presnet paper, we consider quantum variations of the classical random walk 
and refer to such processes as quantum random walks. 

Very recently quantum random walks have been widely investigated by a number of 
groups in connection with the quantum computing, for exmples, [H El IH El El El El UH1 HI] • 
For more general setting including quantum cellular automata, see Meyer This paper is 
an extended version of our previous short letter with no proof (Konno jl]). 

In Ambainis et al. pQ, they gave two general ideas for analyzing quantum random walks. 
One is the path integral approach, the other is the Schrodinger approach. In this paper, we 
take the path integral approach, that is, the probability amplitude of a state for the quantum 
random walk is given as a combinatorial sum over all possible paths leading to that state. 

The quantum random walk considered here is determined by 2 x 2 unitary matrix U stated 
in the next section. The new points of this paper is to introduce 4 matrices, P,Q,R and S 
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given by the unitary matrix U, to obtain a combinatorial expression for the characteristic 
function by using them, and to clarify the dependence of the mth moment and symmetry 
of distribution for the quantum random walk on the unitary matrix U and initial qubit 
(quantum bit) state <p. Furthermore we give a new type of limit theorems for the quantum 
random walk by using our results. Our limit theorem shows that the behavior of quantum 
random walk is remarkably different from that of the classical ramdom walk. As a corollary, 
it reveals whether some simulation results already known are accurate or not. 

The rest of the paper is organized as follows. In Section 2, we introduce a definition of the 
quantum random walk and explain our results. Section 3 gives the characteristic function. 
In Section 4, we present a condition for symmetry of the distribution. Section 5 is devoted 
to a proof of the limit theorem. In Section 6, we consider the Hadamard walk case. 



2 Definition and results. 

The time evolution of the one-dimensional quantum random walk studied here is given by 
the following unitary matrix: 



U = 



a b 
c d 



where a,b,c,d G C and C is the set of the complex numbers. So we have \a\ 2 + |c| 2 = 
\b\ 2 + \d\ 2 — 1, ac + bd — 0, c = — Ab, and d = Aa, where ~z is the complex conjugate of 
zGC and A = det U = ad — be. We note that the unitarity of U gives | A| = 1. 

The quantum random walk is a quantum generalization of the classical random walk in 
one dimension with an additional degree of freedom called the chirality. The chirality takes 
values left and right, and means the direction of the motion of the particle. The evolution of 
the quantum random walk is given by the following way. At each time step, if the particle has 
the left chirality, it moves one unit to the left, and if it has the right chirality, it moves one 
unit to the right. More precisely, we present the left and right chirality states as \L) — *[1, 0] 
and \R) = '[0, 1], where t indicates the transposed operator. So the unitary matrix U acts 
on two chirality states \L) and \R) as 

U\L) = a\L) + c\R), U\R) = b\L) + d\R) . 

The study on the dependence of some important quantities (e.g., characteristic function, 
the mth moment, limit density) on initial qubit state is one of the essential parts, so we 
define the set of initial qubit states as follows: 

$ = {V? = W]GC 2 :M 2 + |/3| 2 = 1}. 

From now on, we will give a precise definition of the quantum random walk X% at time n 
starting from initial qubit state ip G $. First we decompose U = P + Q, where 



P = 



a b 




Q 





c d 



2 



The important point is that P (resp. Q) represents that the particle moves to the left (resp. 
right). We define the (4iV + 2) x (AN + 2) matrix by 



U 



N 






p 












Q 


Q 





p 















Q 





p 


















Q 





p 















Q 





p 


p 












Q 






with 








Then we see that U n becomes also unitary matrix, since P and Q satisfy 



PP* + QQ* = P*P + Q*Q 



1 
1 



PQ* = QP* = Q*P = P*Q 








where * means the adjoint operator. Here an initial qubit state is given by 



On 
On 



E C 



4JV+2 



where Ojv = * [0, . . . , 0] £ C is the zero vector and tp = l [a, The qubit state at time n, 
^^(tp), is determined by 



We write 



¥ n) (<p) 



¥ n \<p) = (u N ) n ¥°\ v ). 



(2.1) 



(n) 



Then ty^\(p) is the two component vector of amplitudes of the particle being at site k and 
at time n with the chirality being left (upper component) and right (lower component). We 
see that (j2.1J) implies 

*i n+1) M = (UN* {n H<p))k = Q^M + P*tiM- (2.2) 

We define the quantum random walk X% with state space {— N, . . . , N} by 

P(XZ = k) = \\9£\<p)\\*. 

By construction P(Xq =0) = 1. We remark P(X% = k) is independent of iV as far as 
n < N. Hence we naturally regard X% as a Z- valued random variable, where Z is the set of 
the integers, and denote by the same symbol X%. 
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In contrast with classical random walks, X% can not be written as X% = Y 1 + • ■ ■ + 
Y n , where Y ± ,Y 2 , . . . are independent and identically distributed random variables. It is 
also noted that the quantum random walk is not a stochastic process. It is a sequence of 
distributions arising from products of the unitary matrix U n- The unitarity of Un ensures 



^P(X* = k) = \\(U N )»V\\ 2 =m\ 2 =\<*\ 2 + 



kez 



for any 1 < n < N and initial state ip = *[0jv, (p, Ojv]. That is, the amplitude always defines 
a probability distribution for the location. For initial state Tp — l [0 N , p,0 N ], we have 



U N (p= [0jv_i,P^,0,Q^,0jv_i], 

V%<p = '[0jv_ 2 , P 2 <p, 0, (PQ + QP)<p, 0, QV, 0jv- 2 ], 

U%Jp = '[On-3, P 3 cp, 0, (P 2 Q + PQP + QP 2 )p, 0, (Q 2 P + QPQ + PQ 2 )y?, 0, Q 3 ip, 0jv_ 3 ]. 



This shows that expansion of C/ n = (P + Q) n for the quantum random walk corresponds to 
that of l n = (p + q) n for the classical random walk. 

We now explain our results briefly. Using an explicit form of P{X% = k) (Lemma 3), we 
obtain the characteristic function of X% (Theorem 4) and the mth moment of it (Corollary 5). 
One of the interesting facts is that, when m is even, the mth moment of X% is independent of 
the initial qubit state G On the other hand, when m is odd, the mth moment depends 
on the initial qubit state. So the standard deviation of X% depends on the initial qubit state 
Our main theorem is as follows: 

THEOREM 1. Suppose abed ^ 0. Then we have 

lim ^ = in law, 

where Z v is a random variable whose distribution has a density f a> p{x)dx such that 



v 1 - \a\ 2 (1 - \ a fix) 2 2 aab(3 + aab(3 

f*A x ) = /i | 2 2 ' mth V/3 = N - \P\ + r-r 2 , 

7r(l — x 2 )^/\a\ 2 — x 2 \o>\ 

for x G (— \a\, \a\), and f a A x ) = ® f or \ x \ — \ a \- Here p = t [a,(3] as before. 

It can be confirmed that f a ,p( x ) satisfies the property of a density function. Indeed, we 
see that f a A x ) > since 1 ± A ajj g|a| > 0, and that 





a 2 


7T 


A/I" 


a 2 


7T 




1. 





f" faA X )dx = ^ ^ [ t- 1 /2 (1 _ t) -l/2 (1 _| a |2 £) -l^ 

y-ioi ;v Jo 

r 

-r(l/2) 2 2 P 1 (l/2,l;l;|o| 2 ) 



Here 2 Fi(a,b; c; z) is the hypergeometric function (see Sect. 5). The last equality comes 
from T(l/2) = v^F and 2 Pi(l/2, 1; 1; \a\ 2 ) = - \a\ 2 . Remark that 



E{Zf) = -(1 - y/T^W) Kp, E((Z*) 2 ) = 1 - JT- 
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Moreover, an easy computation shows that \E((Z ip ) m )\ < 2\a\ m for any m > 1. 

It should be noted that, if \a\ = 1, then b = c = and \d\ = 1. So this case is trivial. 
In fact, Corollary 5 (ii) implies that lim^oo X^/n = W^, in law, where W v is determined 
by PiW^ = —1) = |a| 2 and PiW^ = 1) = |/3| 2 . Theorem 1 suggests the following result on 
symmetry of distribution for the quantum random walk (Theorem 6). Define 

$ s = {<p e $ : P(X% = k) = P(X% = —k) for any n G Z+ and k G Z}, 
$ = {(/? G $ : E(X%) = for any n G Z + } , 

$ ± = = *[«,/?] G $: |a| = |/3| = 1/V2, aabp + aabp = o} , 

where Z + is the set of the positive integers. For (p G $ s , the probability distribution of 
is symmetric for any n G Z + . Using explicit forms of distribution of X% (Lemma 3) and 
E(X%) (Corollary 5 (i) for m — 1 case), we have $ s = $0 = 

3 Characteristic function. 

This section gives a combinatorial expression of the characteristic function of the quantum 
random walk X%. As a corollary, we obtain the mth moment of X%. For fixed / and m, we 
consider 

E(l,m)= P h Q mi P h Q m2 ■■■P l "Q mn . 

ij,mj>0:iiH \-l n =l,mi-\ hm„=m 

It should be noted that for / + m = n and — I + m = k. we see that ^ n) (<^) = E(l,m)(p, since 

(v) = *[^ r il( < ^)' C 2 ) i s a ^ wo component vector of amplitudes of the particle 

being at site k at time n for initial qubit state ip G $ and H(/,m) is the sum of all possible 
paths in the trajectory consisting of I steps left and m steps right with I = (n — k)/2 and 
m — (n + k)/2 (see ()2.2jl ). For example, in the case of P{X\ = —2), we have the following 
expression: 

2(3, 1) = QP 3 + PQP 2 + P 2 QP + P 3 Q. 

Here we find a nice relation: P 2 = aP. By using this, we have S(3, 1) = a 2 QP + aPQP + 
aPQP + a 2 PQ. Moreover, to compute general S(/,m), it is convenient to introduce 

" 

a b 

Then we obtain the following table of products of the matrices P, Q, R and S: 





P 


Q 


R 


S 


p 


aP 


bR 


aR 


bP 


Q 


cS 


dQ 


cQ 


dS 


R 


cP 


dR 


cR 


dP 


S 


aS 


bQ 


aQ 


bS 



R 



c d 




S 
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Table 1 



where PQ = bR, for example. Since P, Q, R and S form an orthonormal basis of the vector 
space of complex 2x2 matrices with respect to the trace inner product (A\B) = tr(A*B), 
S(7,m) has the following form: 

E(l, m) = p n (l, m)P + q n (l, m)Q + r n (l, m)R + s n (l, m)S. 

Next problem is to obtain explicit forms of p n (l,m),q n (l,m),r n (l,'m), and s n (l,m). The 
above example of n = I + m = 4 case, we have 

5(4, 0) = a 3 P, 5(3, 1) = 2abcP + a 2 bR + a 2 cS, 
5(2, 2) = bcdP + abcQ + b(ad + bc)R + c(ad + bc)S, 
5(1, 3) = 2bcdQ + bd 2 R + cd 2 S, 5(2, 2) = d 3 Q. 

So, for example, p 4 (3, 1) = 2abc, 94(3, 1) = 0, r 4 (3, 1) = a 2 b, and s 4 (3, 1) = a 2 c. The following 
holds in general. 

LEMMA 2. We write I Am = min{/,m}. Suppose abed 7^ 0. Then 

(i) for I A m > 1, we have 

H( 1 , m )=^|(-sxr-0(;: 1 1 )[^ + ^-^ + H' 

(ii) for I > 1 and m = 0, we have S(Z, 0) = a l ~ 1 P, 

(Hi) for I = and m > 1, we have 5(0, m) = A m ~ 1 a m ~ 1 Q. 

Proof We first calculate explicit forms of p n (l,m). To begin, we assume / > 2 and m > 1. 
From Table 1, it is sufficient to consider only the following case: 

C(w)^ m = p w ^Q w ^p^ . . .g»2 7 p» 27+ i ; 

where w = (wi) G Z^ 7+1 , with I = J21=o w 2k+i and m = Yll=i w 2k- F° r example, we take 
w± — w 2 — w 3 — 1 and 7 = 1 as PQP. We remark that 27 + 1 is the number of clusters of 
P's and Q's. Next we consider the range of 7. The minimum is 7 = 1, that is, 3 clusters. 
This case is P • • ■ PQ ■ ■ ■ QP ■ ■ ■ P. The maximum is 7 = (l-l)Am. This case includes the 
patterns for example: 

PQPQPQ ■ ■ ■ PQPQPP ■ ■ ■ PP (l-l > m), PQPQPQ ■ ■ ■ PQPQQ ■ ■ ■ QQP (I - 1 < m). 

We introduce a set of sequences with 27 + 1 components: for fixed 7 G [1, (I-l)Am], 

7 7 
W^i, m = {w = (wi) G Z+ 7+1 : ^2w 2k +i = I, ^2w 2k = m}. 

k=0 k=l 

By a standard combinatorial argument, we have 

i^i-C; 1 )^ 1 )- (31) 
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Let w G W yj i !m . Then by using Table 1, we have 

C{w) 1>Um = a^PcT^Qa^P ■ ■ ■ (T^Qa^+^P 
= a l -^ +1) d m ~' f (PQyP 

= a l -^ +l) d m ^b' 1 c<P. (3.2) 
Combining (|3.1|) and (|3.2j) . we obtain 

(i-l)Am 

p n (l,m)P = ^2 ^ C { w )i^,m 
7=1 ™ew 7iiim 

(l-l)Am 



7=1 

When I > 1 and m = 0, it is easy to see that p n (l,0)P = P l = a l ~ 1 P. Furthermore, when 
/ = 1, m > 1 and I — 0, m > 0, it is clear that p n (/, m) = 0. 

As in the case of p n (l,m), we compute q n (l,m),r n (l,m), and s n (l,m) by considering the 
patterns Q Wl P W2 Q W3 ■ ■ ■ p^q^+i^ p Wl Q W2 p WA . . . qw 2 ^ and Qmp^gm . . . respec _ 

tively. Then we obtain 

£|A(m-l) ^ m -i^ a/ _ 767c7rfm _ (7+1) for / > 1; m > 2; 

g n (7,m)=<( rf™- 1 for/ = 0,m>l, 

for m = 1, Z > 1 and m = 0, Z > 0, 

E'ti ( 7 _i) Ha^ 1 ^ for Z, m > 1, 



jj=i V7— 1/V7— 1/~ "~ ~ — , 

for / A m = 0, 

E7I1 ( 7 li) (^o'-Tft-r-Vd^ for Z, m > 1, 

for Z A m = 0. 



r n (l,m) = 
s n (l,m) = 

For Z A m > 1, the above explicit forms of p(Z, m), g n (Z, m), r n (Z, m), and s n (Z, m) imply 



07 07 c 



From c = — Ab and d = Aa, the proof of Lemma 2 (i) is complete. Furthermore, parts (ii) 
and (iii) are easily shown, so we will omit the proofs of them. □ 

The distribution of X% can be derived from Lemma 121 by direct computation. Let [x] 
denote the maximal integer smaller than or equal to x. Let 



^7,5,n,fc 



k - 1\ fk - 1\ fn - k - 1\ (n - k - 1 
7-lJU-lM 7 -l M 5-1 
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LEMMA 3. For k = 1,2,..., [n/2], we have 
P(X% = n-2k) 



k k 

■i 21 """ EE 

7=1 5=1 



H2 \ 7+<5 



i a r/ V 7^ 
+{A;>| 2 + (n-A;) 2 |a| 2 -(7 + 5)A;}|/5| 2 



{k 2 \a\ 2 + (n - fc) 2 |6| 2 - ( 7 + - fc)}|a| 



+ 



|6P 



{(n - % - k5 + n(2k - n)\b\ 2 }aabp 



P(X* = -{n-2k)) 

k k 

= h 2( " _1) EE 

7=1 8=1 



+{-^7 + (n - + n(2k - n)|6| 2 }aafe/3 + 7<5 

| 6 |2\7+* 



^"y,5,n,k 
7(5 



{A; 2 |6| 2 + (n-A;) 2 |a| 2 -(7 + 5)A;}|« 



+{fc 2 |a| 2 + (n - fc) 2 |6| 2 - (7 + <T)(n - A;)}|/5| 2 

+ T7T7 1"{^7 - (n - - n(2k - n)|&| 2 }aa<&/3 
|fo| 2 L 

+{-(n - fc)7 + fc5 - n(2A; - n)|6| 2 }aa6/? + 7<5 

P(X^ = n) = H^^Hl&rM 2 + \a\ 2 \P\ 2 - (aabp + aabp)}, 
P(X% = -n) = \a\ 2{n - l) {\a\ 2 \a\ 2 + \b\ 2 \(3\ 2 + {aabp + aab(3)} . 

By using Lemma^l we obtain a combinatorial expression for the characteristic function of 
X% as follows. This result will be used in order to obtain a limit theorem of X%. Let /i Q / g = 

(H 2 - \b\ 2 ) (H 2 - \(3\ 2 ) + 2(aab]3 + aab(3) and u^ s ,n,k = (n- k) 2 + k 2 - 71(7 + 5) + 2 7 <v'H 2 . 

THEOREM 4. (ij Suppose abed + 0. T/ien we have 



E(e 



| a |2(«-l) 



EEE 

fc=i 7=1 5=1 



cos(n£) — z/i aj/ 3 sin(n^) 

| 6 |2X7+5 



7<5 



^7,5,n,fe cos((n - 2k)£) 



-(n - 2A;)<{ n + -^tttt(|«| 2 - |/?| 2 - V<x,p) \isax{{n - 2k)£) 



n n 
2 2 



-EE 

7=1 <5=1 



|6| 



2|6| 2 

2 \ 7+<5 



^7, (5,n, n/2 

2 7 5 



^7, (5, n, n/2 



where I(x) = 1 (resp. =0) if x = (resp. x 7^ 0). 
(wj Lei 6 = 0. JTien we /lave 



£(e^ x " ) = cos«) + z(|/3| 2 - |a| 2 ) sin«). 
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(Hi) Let a = 0. Then we have 

E( e ^ x ") = { cos ^ + i (\ a \ 2 ~ l/?| 2 ) sin £ l f n is °dd, 
\ 1 if n is even. 

We should remark that the above expression of the characteristic function in part (i) is not 
uniquely determined. From this theorem, we have the mth moment of X^ in the standard 
fashion. The following result can be used in order to study symmetry of distribution of X%. 

COROLLARY 5. (i) Suppose abed ^ 0. When m is odd, we have 



\n— 11 



<P\m\ _| |2(n-l) 



k=l 7=1 ,5=1 



- J — — / \b\ 2 V +5 {n-2kT^K lAn ,k 



78 



rn-il 



When m is even, we have 

2 J k k , i, |2 \ 7+5 
fe=l 7 =1 5=1 



(ii) Let 6 = 0. Then we have 

e{{KT) 



n m (\(3\ 2 — | or | 2 ) if m is odd, 
n m if m is even. 



(Hi) Let a = 0. Then we have 



\a\ 2 — \/3\ 2 if n and m are odd, 
E((X%) m ) = ^ 1 if n is odd and m is even, 

if n is even. 

For any case, when m is even, E((X%) m ) is independent of initial qubit state </?. Therefore 
a parity law of the mth moment can be derived from the above result. 



4 Symmetry of distribution. 

In this section, we give a necessary and sufficient condition for the symmetry of the distri- 
bution of X%. 

THEOREM 6. Let $ s ,$ 0; an d $± be as in Section 2. Suppose abed ^ 0. Then we have 
$ s = $ = $ ± . 

This is a generalization of the result given by [3] for the Hadamard walk introduced in 
Section 6. Nayak and Vishwanath jHj discussed the symmetry of distribution and showed 
that *[1/V2, ±i/V2] e $ s for the Hadamard walk. 
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Proof, (i) $ s C $o- This is obvious by definition. 

(ii) $ C By Corollary 5 (i) with m = 1, we see that E(Xf) = E{X$) = if and only 
if A*o,/3 = 0- Then this implies that for n > 3, Corollary 5 (i) with m = 1 can be rewritten as 

[2=11 



eVv^ l»l * HI"! ~ i^m ; \^ \^ \^ 



2 |6|2 | a |2/ g 

I I fc=l 7 =1 5=1 V 1 1 ' ' 

Therefore E(X%) = (n > 3) gives |a| = |/3|. Combining |a| = \/3\ with /i aj/ g = 0, we have 
the desired result. 

(iii) $j_ C $ s . We assume that \a\ = |/3| = l/v2 and aa6/3 + aab(3 = 0. By using these and 
Lemma El we see that for k = 1, 2, ... , [n/2], 

I |2(n— 1) k k , |/ 1 |2\7+^ 

W = n - 2*) = P« = -(n - 2*)) = M— £ £ ( "Ai ) 

/ 7 =i a=i V l a l / 7f 

and P(X% = n) = P(X% = —n) = jap 71-1 ) \a\ 2 . So the desired conclusion is obtained. □ 



5 Proof of Theorem 1. 

Let P^ ,fl (x) denote the Jacobi polynomial. Then it is well known that P^'^ix) is orthogonal 
on [—1, 1] with respect to (1 — x) u (l + x) M with is, fi > — 1, and that the following relation 
holds: 

P ^ {X) = r(n + l)r(t+l) 2iri( ~ W ' n + ^ + ^ + 1 ^ + 1 ^ 1 ~ z)/2), (5.1) 

where T(z) is the gamma function and 2F±(a, b; c; z) is the hypergeometric function: 

^ T(a + n)T(b + n) T(c) z« 

2 F 1 {a, b c;z)= > — — - • — . 

/ — J 1 (a) 1 (6) 1 (c + n) n! 



n=0 



Let p n ,A;,i = Pfc-i 2fc (2 1 a | 2 — 1) for i = 0,1. Then we see that 

1 fk - 1\ /ra - Jfe - 
7 VT - V V 7-1 



7=1 



-2(fc-i) 2Fi( _ (A . _ 1), n _ A; + 1; 2; 1 - |o| 



2\ 



The first equality is given by the definition of the hypergeometric function (see p. 35 in [TT]). 
The second equality comes from the relation: 2 Pi(o, 6; c; = (1— ;z)~ a 2Pi(a, c— 6; c; z/(z— 1)). 
The last equality follows from (|5.1j) . In a similar way, we have 



k / |6| 2 \ 7_1 A-l\ fn-k-1 



£ 

7=1 



7 ~ 1/ \ 7 — 1 



-2(fc-l)„ 

Pn,k,0- 



By using the above relations and Theorem 4, we obtain the following asymptotics of char- 
acteristic function E(e l ^ Xn ^ n ): 
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LEMMA 7. Ifn^oo with k/n = x G (-(1 - \a\)/2, (1 + \a\)/2), then 

E{e* x -' n ) ~ ^ |a| 2 "- 4fe - 2 |6| 4 

2 2 1 

- - Pn,kfiPn,k,l + Pn,fe,0 j «)s((l - 2x)£) 



k=l 

2x 2 -2x + l 



X 2 



- ( ~ J {— ^,1 + p n , fc ,oPn,fc,i]«sm((l -2x)0 

where fin) ~ g(n) means f{n)/g(n) — > 1 (n — > oo). 

Next we use an asymptotic result on the Jacobi polynomial p^+ an ^+ hn ^ x ) derived by 
Chen and Ismail [2|. By using (2.16) in their paper with a -> or 1, a - > 0, (3 = b — > 
(1 — 2x)/x,x — > 2\a\ 2 — 1 and A — > 4(1 — |a| 2 ){(2x — l) 2 — \a\ 2 }/x 2 , we have the following 
lemma. It should be noted that there are some minor errors in (2.16) in that paper, for 

example, a/ (—A) — > \J (—A) 

LEMMA 8. Ifn^oo with k/n = x G (-(1 - \a\)/2, (1 + \a\)/2), then 

2 1 „ 1 2A; — ri 

p nA o . cos(An + B), 

V imy—K. 

p nA i - . / - j-^ cos(An + B + 6), 



where A = (1 — |a| 2 ){(2x — l) 2 — |a| 2 }, A and I? are some constants (which are independent 
of n), and 9 G [0, tt/2] is determined by cos 9 = *JJT— |a| 2 )/4x(l — x). 

Proof of Theorem 1. From the Riemann-Lebesgue lemma and Lemmas [7| and |H1 we see 
that 

limE(e^) = ^—Mi [~ 2 cos((l - 2x)j) - i\ aS {\ - 2x) sin((l - 2x)Q ^ 
n ^°° ' 7T yidai x (\ _ x )y/(\a\ 2 - l)(4a; 2 - Ax + 1 - \a\ 2 ) 



a/1 — |a| 2 cos(x£) — i\ a ,p x sin(x£) 



7T 



|a| (1 - X 2 ) A/|a| 2 - X 2 



\a\ 



^1 - \a\ 2 (1- X a ^x) ^ x ^ 
\a\ 7r(l - x 2 )\/|a| 2 - X 2 
Hence X%/n converges weakly to the limit Z v . □ 



6 Hadamard walk case. 



In this section, we focus on the Hadamard walk, which has been extensively investigated in 
the study of quantum random walks. The unitary matrix U of the Hadamard walk is defined 
by the following Hadamard gate (see Nielsen and Chuang [TUj): 



U 



1 

V2 
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The dynamics of this walk corresponds to that of the symmetric random walk in the classical 
case. However the symmetry of the walk depends heavily on initial qubit state, see [S]. 

For example, in the case of the Hadamard walk with initial qubit state (p = * [1/ \/2, i/\/2] 
(symmetric case), direct computation gives 

P{Xt = -4) = P{Xt = 4) = 1/16, 

P(Xf = -2) = P(X% = 2) = 6/16, P(X$ = 0) = 2/16. 

In contrast with the above result, as for the classical symmetric random walk Y° starting 
from the origin, we see that 

P(Y 4 ° = -4) = P(Y° = 4) = 1/16, 

P(F 4 ° = -2) = P{Y° = 2) = 4/16, P(Y° = 0) = 6/16. 

In fact, quantum random walks behave quite differently from classical random walks. For 
the classical walk, the probability distribution is a binomial distribution. On the other hand, 
the probability distribution in the quantum random walk has a complicated and oscillatory 
form. 

Now we compare our analytical result (Theorem 1) with the numerical ones given by 
Mackay et al. jO], Travaglione and Milburn ^3] for the Hadamard walk. In this case, 
Theorem 1 implies that for any initial qubit state (p — t [a, 

hm P(a < X*/n < b) = | n{1 _ X 2 )VT ^ VvaW*) dx ' 

where l/ u>v )(x) is the indicator function, that is, l( u> v){%) = 1, if ic G (u, v), — 0, if x ^ (u,v). 
For the classical symmetric random walk Y° starting from the origin, the de Moivre-Laplace 
theorem shows 

rb e -x 2 /2 

lim P(a < Y°/^ <b)= ^=dx. 

' J a V27T 

If we take (p = *[l/\/2, i/y/2) (symmetric case), then we have the following quantum version 
of the de Moivre-Laplace theorem: 

nb l 

lim P(a < X'f/n < b) = / 1, lA /o, / /nAx) dx. 

™ V " n ' ~ ' J a 7T(1 - X 2 ) VI - 2X 2 (-W2,W2)V J 

So there is a remarkable difference between the quantum random walk X% and the classical 
one Y° even in a symmetric case. Noting that E(X%) = (n > 0) for any y9 G $ s , we have 



lim sd(X%)/n = J (2 - V2)/2 = 0.54119 . . . , 

where sd(X) is the standard deviation of X. This rigorous result reveals that numerical 
simulation result 3/5 = 0.6 given by [13] is not so accurate. 
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As in a similar way, if we take tp = *[0, e ld ] with 9 e [0, 2tc) (asymmetric case), then we 

see 

f b 1 

lim P(a < X^ln < b) = / 1, , , r~ , , /^.(x) dx. 

™ V ~ nl ~ 1 J a tt(1 - x)y/l - 2x 2 (-V^V^V ; 

So we have 

lim E(X%)/n = (2 - = 0.29289 . . . , lim sd(X^)/n = J {y/2 - l)/2 = 0.45508 . . . . 

n^oo n— >oo " 

When <yj = *[0, 1] (0 = 0), Nayak and Vishwanath 9J and Ambainis et a/. P gave a 
similar result, but both papers did not treat weak convergence. The former paper took the 
Schrodinger approach, and the latter paper took two approaches, that is, the Schrodinger 
approach and the path integral approach. However both their results come mainly from the 
Schrodinger approach by using a Fourier analysis. The details on the derivation based on 
the path integral approach in [1] are not so clear compared with this paper. 

In another asymmetric case (p = t [e , 0] with 9 e [0, 2ir), a similar argument implies 

f b 1 

lim P(a < X'f/n < b) = / 1, , , ^, , ^(x) dx. 

™ V " n ' ~ ' Ja 7T(1+X)v / 1-2X 2 (-W2,1A/2)V ) 

The symmetry of distribution gives the following same result as in the previous case (p = 

t [0,e l9 }. So the standard deviation of the limit distribution Z v is given by yj (y/2 — l)/2 = 
0.45508 .... Simulation result 0.4544 ± 0.0012 in [Oj (their case is 9 = 0) is consistent with 
our rigorous result. 

Acknowledgment. The author would like to thank the referee for the careful reading and 
useful suggestions which improve the paper. 
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